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Abstract. We perform a study of the ultra high energy neutrino detection 
performances of a km 3 Neutrino Telescope sitting at the three proposed sites for 
ANTARES, NEMO and NESTOR in the Mediterranean sea. Wc focus on the effect of the 
underwater surface profile on the total amount of yearly expected r and /i crossing 
the fiducial volume in the limit of full detection efficiency and energy resolution. Wc 
also emphasize the possible enhancement of matter effect by a suitable choice of the 
geometry of the Telescope. 
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1. Introduction 

Neutrinos are one of the main components of the cosmic radiation in the ultra-high 
energy (UHE) regime. Although their fluxes are uncertain and depend on the production 
mechanism, their detection can provide information on the sources and origin of the UHE 
cosmic rays. For example, UHE neutrinos can be produced via 7r-photoproduction by 
strongly accelerated hadrons in presence of a background electromagnetic field. This 
scenario is expected to occur in extreme astrophysical environments like the jets of active 
galactic nuclei, radio galaxies and gamma ray burst sources as well as in the propagation 
of UHE nucleons scattering off the cosmic background radiation (known as cosmogenic 
neutrinos [H [2]). 

From the experimental point of view, after the first pioneering and successful 
achievements, neutrino astronomy in the high energy regime [31 HI El El E] is a 
rapidly developing field, with a new generation of neutrino telescopes on the way. 
A benchmark result was obtained by the DUMAND [5] collaboration, followed by the 
successful deployments of NT-200 at Lake Baikal [9] and AMANDA [10] at the South 
Pole, which have shown the feasibility of large optical Cherenkov neutrino telescopes 
(NT) in open media like sea- or lake-water and glacial ice. These experiments observed 
atmospheric neutrinos [TTj and set bounds on their extraterrestrial flux [121 (131 EH] which 
are much more constraining than the corresponding bounds obtained by underground 
neutrino detectors [15]. These interesting results and the perspective to perform 
astronomical studies using UHE neutrinos stimulated several proposals and R&D 
projects for neutrino telescopes in the deep water of the Mediterranean sea, namely 
ANTARES P], NESTOR [HI EH OS] and NEMO [20], which in the future could lead to the 
construction of a km 3 telescope as pursued by the KM3NeT project [2H[22]. Actually, the 
ANTARES collaboration is in a more advanced phase, with a telescope with an area of ~ 0.1 
km 2 already under construction [23]. A further project is IceCube, a cubic-kilometer 
under-ice neutrino detector [2l"[ [25l [26] currently being deployed in a location near 
the geographic South Pole in Antarctica. IceCube applies and improves the successful 
technique of AMANDA to a larger volume. 

Until now, the possibility to perform astronomy with neutrinos has been seriously 
limited by the presence of the heavy atmospheric background in the energy range 
presently explored. The start of //-astronomy is eagerly waiting the completion of 
the new km 3 projects. In fact, according to theoretical expectations, the km 3 is the 
minimum detector size required to detect with reasonable chances of success point-like 
sources at the TeV scale and, more relevant for this paper, to explore energies above 
about 100 TeV, where extraterrestrial diffuse fluxes should start to dominate over the 
steeper atmospheric spectrum. In the following we shall focus the attention mainly on 
this energy range, although many of our results are valid also at lower energies. 

Although NT's were originally thought as detectors, their capability as v T 
detectors has become a hot topic [271 EH ESI ffl EH [32l [33l [341 [35] , in view of the fact 
that neutrino oscillations lead to nearly equal astrophysical fluxes for the three neutrino 
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flavorqj. Despite the different behavior of the produced tau leptons with respect to 
muons in terms of energy loss and decay length, both and v T detection are sensitive 
to the matter distribution near the NT site. Thus, a computation of the event detection 
rate of a km 3 telescope requires a careful analysis of the surroundings of the proposed 
site. The importance of the elevation profile of the Earth surface around the detector 
was already found of some relevance in Ref. [H], where some of the present authors 
calculated the aperture of the Pierre Auger Observatory [121 H3] for Earth-skimming 
UHE z/ r 's. Indeed, air shower experiments can be used as NT's at energies > 10 18 eV, a 
topic recently reviewed in [H]. In particular, the possibility of detection of the r leptons 
produced by Earth-skimming UHE z/ T 's has been analyzed in a series of papers [H], [^5] - 
[55] . In Ref. [H] the use of a Digital Elevation Map (DEM) of the geographical area of 
the experiment proved useful to characterize peculiar matter effects in Earth-skimming 
events. 

The aim of this paper is to estimate the effective aperture for v T and detection of 
a km 3 NT in the Mediterranean sea placed at any of the three locations proposed by the 
ANTARES, NEMO and NESTOR collaborations. We do not consider any detail related to the 
experimental setup nor the detector response. In particular, we assume full detection 
efficiency via Cherenkov radiation for muons and taus crossing the NT fiducial volume. 
This means that we consider a lepton as detected if it crosses in any point the surface 
delimiting our fiducial volume, without e.g. taking into account further requirements 
or cuts needed for a good directional or energy reconstruction. These would depend on 
several parameters, like the spacing of the strings, the distribution of photomultipliers 
along the string, etc. which are characteristics of the apparatus and thus beyond the aim 
of the present analysis. We rather compare the site characteristics by using the DEM 
of the different areas. We shall therefore characterize and quantify the importance of 
"matter effects" for the three sites, and focus on the role played by the geometry of 
the experiment in enhancing the effect. These considerations may provide an important 
ingredient in shaping the final design of a km 3 Mediterranean NT. 

A detailed DEM of the under-water Earth surface is available from the Global 
Relief Data survey (ETOP02) [56], a grid of altimetry measurements with a vertical 
resolution of 1 m averaged over cells of 2 minutes of latitude and longitude. In Figures 
[Q [2] and [3] we show the 3D maps of the areas around the three NT sites. The black curve 
represents the coast line, whereas the red spot stands for the location of the apparatus. 
By following the same approach developed in [H], we use this DEM to produce a realistic 
and statistically significant sample of v T jr and tracks crossing the fiducial volume 
of the NT that are then used to evaluate the effective aperture of each detector. 

We note that when the events are reconstructed only in terms of the energy loss 
along the track, the UHE taus can not be distinguished from less energetic muons. This 
implies that the reconstruction analyses of UHE and v T events are highly entangled 
issues. We shall consider both of them, although for the sake of clarity we shall first 

% This statement may not hold for exotic neutrino models [361 EI] or f° r peculiar astrophysical sources 
[381 [391 [40]. 
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Figure 1. The surface profile of the area near the ANTARES site (red spot) at 42° 30 
N, 07° 00 E. The black curve represents the coast line. The sea plateau depth in the 
simulation is assumed to be 2685 m. The effective volume starts at an height of 100 
m from the seabed, to account for the spacing of the first photomultipliers as foresee 
by the current designs. The km 3 detector is oriented along the E-W/S-N directions. 



focus on v T detection. Of course, when considering shower events or more in general 
contained events — where the charged lepton production happens inside the instrumented 
volume — including events with peculiar topologies like "lollipop" or "double bang" , there 
are realistic chances of flavor-tagging in the detector. However, in the UHE range above 
~ 10 7 GeV these kind of events are subdominant with respect to the bulk of tau track- 
events. Further details on the flavor discrimination possibilities are discussed in Section 



The structure of the paper is as follows. In Section [2] we introduce the formalism 
and definitions used in the analysis and define the aperture for a NT. Our results for 
v T induced events are reported and discussed in Section [3] for various incoming neutrino 
fluxes, while v^jn events are described in Section HI Finally, we report our conclusions 
in Section [51 

2. The effective aperture of a NT 

We define the km 3 NT fiducial volume as that bounded by the six lateral surfaces S a 
(the subindex a=D, U, S, N, W, and E labels each surface through its orientation: 
Down, Up, South, North, West, and East), and indicate with Vt a = (9 a ,(fi a ) the generic 
direction of a track entering the surface S a . The scheme of the NT fiducial volume 
and two examples of incoming tracks are shown in Fig. HI We introduce all relevant 
quantities with reference to v T events, the case of being completely analogous. 

Let d$ u /(dE u dQ a ) be the differential flux of UHE v T + v T . The number per unit 



HI 
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Figure 2. The surface profile of the area near the NEMO site (red spot) at 36° 21 N, 
16° 10 E. The black curve represents the coast line. The sea plateau depth used in 
the simulation is 3424 m. The effective volume starts at an height of 100 m from the 
seabed, to account for the spacing of the first photomultipliers as foresee by the current 
designs. The km 3 detector is oriented along the E-W/S-N directions. 



time of r leptons emerging from the Earth surface and entering the NT through S a with 
energy E T is given by 

( dN A [aC> f A Q ( ATP d *»(Eu, »a) 

x I dE T cos{e a )k T a {E u ,E T ;f a ,Q a ) . (1) 



This equation is the same as in [?T], but for full duty cycle and detection efficiency 
The kernel k T a [E v , E T ; f a , fl a ) is the probability that an incoming v T crossing the Earth, 
with energy E v and direction Q a , produces a r-lepton which enters the NT fiducial 
volume through the lateral surface dS a at the position r a with energy E T (see Fig. 0] for 
the angle definition). If we split the possible events between those with track intersecting 
the rock and the ones only crossing water, the kernel k T a {E v , E T ; f a , fl a ) is given by the 
sum of these two mutually exclusive contributions, 

kl(E„ , E T ; fa, fi„) = k T /(E u , E T ; f a , fi a ) + k T a w (E u , E T ; f a , fi n ) . (2) 

Let us focus on the rock events contributing to k% r (E v , E T ;r a ,Q a ). These can be 
classified according to their production mechanism as follows: 

1) events in which the v T interacts producing a r in the rock (rl); 

2) events in which the v T interacts producing a r in water, on the way to the NT (r2); 

3) events in which the v T interacts producing a r inside the NT fiducial volume (r3). 
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Figure 3. The surface profile of the area near the NESTOR site (red spot) at 36° 21 
N, 21° 21 E. The black curve represents the coast line. The sea plateau depth in the 
simulation is assumed to be 4166 m. The effective volume starts at an height of 100 
m from the seabed, to account for the spacing of the first photomultipliers as foresee 
by the current designs. The km 3 detector is oriented along the E-W/S-N directions. 




Therefore one has 

k T /{E v , E T ; f a , Sl a ) = k T /\E v , E T ; f a , fi a ) + k T / 2 (E y , E T ; f a , ft a ) 

+ k:' r3 (E u ,E T ;f a ,Q a ) . (3) 

Although here, for the sake of brevity, we only discuss in details the events occurring 
in rock (rl), the analysis of those of type (r2) and (r3) is completely analogous and 
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straightforward. Of course, all contributions (rl)-(r3) have been added to compute the 
event rate. 

As already shown in details in [^Tl 153] a (rl)-event corresponds to the simultaneous 
fulfillment of the following conditions: 

1) A v T with energy E v travels over a distance z through the Earth before interacting. 
The corresponding probability P\ is given by 



bd^l' (4) 



where 



where Na the Avogadro number. See [HI [53] for notations as well as a detailed 
discussion of the neutrino-nucleon cross section, (7qq(E v ). In the present formalism 
the effect of the Earth density profile is approximated using, track by track, the 
averaged g r along the chord subtended by that track. The calculations are made 
with the parametrization of the Earth density profile of Ref . [57] . Note, however, 
that for almost horizontal events particles travel in the terrestrial crust only, and 
thus the superficial value for the Earth density g r ~ 2.65 g/cm 3 [51] would be an 
accurate approximation. Some differences could appear for low energy particles 
deeply crossing the Earth. We checked that, using the constant value of the crust 
density for all the Earth density profile the changes in the tau aperture are generally 
less than 10%, while the effect is of the order of 16% for muons with E < 10 5 GeV, 
i.e. below the energy range for which the cosmic flux is expected to dominate 
over the atmospheric flux. The inclusion of the Earth density profile also affect at 
the 10% level the distributions of both rs and fis incoming zenith angles: due to 
increased screening along the nadir direction (vertical up-going), the distributions 
slightly shrink along the horizontal direction. 

2) The neutrino produces a r in the interval z, z + dz, the probability of such an event 
being 

P 2 dz = xv dz . (6) 

*CC \ E v) 

We do not consider here the event corresponding to the scattering of a v T via 
neutral current in the Earth followed by conversion via charged current, which 
amounts to a small distortion of the incoming neutrino flux, the latter being yet 
unknown [59]. Of course, this sub- leading effect should be added when trying to 
reconstruct the flux from experimental data. Also, we consider the charged lepton 
track as collinear with the parent neutrino direction, which is highly accurate given 
the huge relativistic boosting factors involved. 

3) The produced r emerges from the Earth rock with an energy E' T . This happens 
with a probability 

f II 1 \ I 
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x 5 (E' T - E Q T (E U ) e -/W(*-*)) , (7) 

where m T = 1.78 GeV, r T ~ 3.4 x 10~ 13 s is the r mean lifetime and E® is the r 
energy at production, whereas the parameter j3 T = 0.71 x 10~ 6 cm 2 g _1 weights the 
leading term in the r differential energy loss in rock [531 160] . 
dE 

-r^ = -((3 T + lT E T )E T g r . (8) 
dz 

The contribution of 7 r can be neglected as it only affects extremely energetic r's 
which, differently from the case of the Pierre Auger Observatory, are not relevant 
for NT's. The quantity z r (r a , Q a ) represents the total length in rock for a given 
track entering the lateral surface S a of the fiducial volume at the point r a and with 
direction fl a . 

4) Finally, the r lepton emerging from the Earth rock propagates in water and enters 
the NT fiducial volume through the lateral surface S a at the point f a with energy 
E T . The corresponding survival probability is 

0) 

where g w stands for the water density and z w (r a , Q a ) represents the total length 
in water before arriving to the fiducial volume for a given track entering the lateral 
surface S a at the point f a and with direction Q a . 

Collecting together the different probabilities in Eqs. flU), (JSD, dZl) and Q, we have 

f zr rE° T (E v ) 

kY\E v , E T ■ f Q , n a ) = dz dE' T Pi P 2 P 3 P 4 . (10) 

Jo Jo 

Similar results can be obtained for the (r2)- and (r3)-events as well as for those we 
defined as water-like. 

For an isotropic flux we can rewrite Eq. (0Q), summing over all the surfaces, as 

dt J 47r dE u 

which defines the total aperture A T ^ r,w \Ey), with "r" and 'W denoting the rock and 
water kind of events, respectively. The contribution of each surface to the total aperture 
reads 

A£ r > w \E v ) = JdE T J dQ a J dS a cos (9 a ) k T / r ' w \E u , E T ; f a , Q a ) . (12) 
3. The event rate for v T interactions 



We show in Fig. [5] the apertures A T ^ r,w ' > for the NEMO site together with the corresponding 
quantity for the Pierre Auger Observatory Fluorescence Detector (FD) calculated in |41j . 
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Figure 5. The effective apertures A t{ - t \E v ) (solid line) and A t{ - w) (E u ) (dashed line) 
defined in Eq. Ijlip versus tau neutrino energy for NEMO. The dotted line corresponds 
to the same quantity for the Auger Fluorescence Detector for Earth-skimming v r as 
in gl]. 




Figure 6. The effective apertures A T a (E v ) of Eq. \12\ versus tau neutrino energy 
for (left) rock events and (right) water events for the NEMO site. 



Note that the Auger case is only for Earth-skimming r's, since down-going neutrino 
induced events can be disentangled from ordinary cosmic rays only for very inclined 
showers. Interestingly, the NEMO- water and Auger-FD apertures almost match at the 
FD threshold of 10 18 eV, so that using both detectors results into a wide energy range 
of sensitivity to v T fluxes. 
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Figure 7. Contour plot in the plane zenith angle-r energy for the NEMO site and for 
rock (red full lines) and water (black dashed lines) events. In both cases the contours 
enclose 68, 95 and 99 % of the total number of events calculated assuming a GZK- 
WB flux ([HI], see also [SHj). cos 8 = 1,0,-1 correspond respectively to down-going, 
earth- skimming, and up-going events. 

We show in Fig. [6] the high energy behavior for each surface contributing to the 
effective aperture. For rock events there is a clear W-E asymmetry, easily understood 
in terms of matter effects related to the particular morphology of the NEMO site (see Fig. 
[2]). A much smaller S-N asymmetry is also present. In other words, the asymmetries in 
the number of rock-events reflect the asymmetries in the morphology of the site. 

For neutrino energies larger than 10 7 GeV the main contribution to the aperture 
comes from the lateral surfaces, i.e. from r leptons emerging from the rock 
far from the NT basis and crossing the fiducial volume almost horizontally Instead, 
the upper surface contribution is negligible due to the very small fraction of events 
crossing the rock and entering the detector from above. The decreasing contribution of 
the bottom face to rock events is due to the Earth shadowing effect. 

For water events the contribution to the aperture from all surfaces is comparable 
(except for the lower one which has no events), the upper one providing a slightly larger 
contribution as the energy decreases. Indeed, events which would cross the lateral 
surfaces should travel over a longer path in water and this becomes more unlikely at 
lower energies due to the shorter r decay length. The matter effect in the case of water 
events is less pronounced and anticorrelated with the asymmetries in the morphology 
of the site resulting effectively in a small (percent) screening effect. 

In Fig. [7] we report, for both the rock and water cases and for the NEMO site, the 
contours enclosing 68, 95 and 99 % of the total event rate, as they appear in the 
plane E T -9 plane, where 9 is the arrival direction zenith angle. These results were 
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Figure 8. A comparison of the effective apertures A T ^ r ' w ^(E v ) for the three NT 
sites. We plot the ratios [A T » r ' t0 » (NESTOR) - A T ^ r ' w ^(mm)]/A T ^ w ^(mm) and 
[^((r,™)) (ANTARES) - A T (( r ' tu »(NEM0)]/A T » r < tu »(NEM0) versus the neutrino energy. 



obtained assuming a Waxmann-Bahcall-like neutrino flux (GZK-WB) [6T] (see also [JT] 
and references therein). As the energy increases, the arrival directions of rock events are 
almost restricted to the horizontal (Earth-skimming), while at lower energies the earth- 
screening effect is less pronounced and this explains the broader angular distribution. 
The situation is different for water events for which the angular distribution is broad 
at all energies, a purely geometrical effect due to the fact that (down-going) water 
events are not screened in few kilometers of water. The same geometrical considerations 
explain the ratio of water to rock event rate of O(10) (see Table [T|) which is simply 
related to the ratio of down-going to Earth-skimming solid angles. This is the same 
kind of behavior expected in the Auger detector although, as we already mentioned, the 
down-going events in this case are hardly distinguishable from the background of proton- 
induced showers so that only Earth-skimming or almost horizontal showers can be used 
to identify unambiguously the neutrino-induced events. It is also worth commenting the 
expected r energy distribution shown in Fig. [7J All events correspond to a relatively 
narrow energy window, from 10 6 GeV up to 10 10 GeV, where the lower cut-off arises 
from the shorter r decay length at low energy. 

In Fig. [8] we compare the detection performances of a km 3 NT placed at one of the 
three sites in the Mediterranean sea. The NESTOR site shows the highest values of the 
r-aperture for both rock and water, due to its larger depth and the particular matter 
distribution of the surrounding area, while the lowest rates are obtained for ANTARES. 
The aperture in the three sites can be quite different at high energy, but in order to get 
the expected number of UHE events per year, one has to convolve the aperture with a 
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0002/0 2133 


0002 /0 2543 


s 


0.0185/0.1602 


0.0256/0.1773 


0.0240/0.2011 


N 


0.0241/0.1540 


0.0229/0.1823 


0.0321/0.1924 


W 


0.0212/0.1584 


0.0335/0.1691 


0.0265/0.2002 


E 


0.0206/0.1589 


0.0190/0.1875 


0.0348/0.1907 


Total 


0.090/0.799 


0.107/0.929 


0.123/1.039 



Table 1. Estimated rate per year of rock/water r events at the three km 3 NT sites 
for a GZK-WB flux ([STJ[SS]). The contribution of each detector surface to the total 
number of events is also reported. 

neutrino flux which typically drops rapidly with the energy. Although the percentage 
value of the matter effects remains unchanged, in this very low statistics regime they 
can be hardly distinguished; still, they can be enhanced by an appropriate choice of the 
detector shape as we discuss in the following. 

Knowing the aperture of the NT at each site, we can compute the expected r event 
rate, once a neutrino flux is specified. In Table [T]these rates are shown assuming a GZK- 
WB flux ( (6TJ [55]). The effect due to the local matter distribution is responsible for the 
N-S, W-E and NE-SW asymmetries for the ANTARES, NEMO and NESTOR sites, respectively, 
as expected from the matter profiles shown in Figs. [1], [2] and [31 These matter effects, 
for the specific UHE flux considered (GZK-WB), correspond to an enhancement of rock 
events which goes from 20 to 50% for the three sites, respectively, and a screening factor 
for water events from 3 to 10%. The largest relative difference among lateral surfaces 
is in the case of W/E for NEMO, where the huge wall to the west of the site (see Fig. [2]) 
improves the rate by about 75%, almost a factor 2! Notice also that the water events 
from the U surface are basically proportional to the depth. 

It is important to emphasize that the impact of matter effects on the rates depends 
critically upon the energy spectrum of the UHE neutrino flux. For more energetic 
neutrino fluxes the enhancement factor is expected to be more significant (see the energy 
behavior of A T J- r \E v ) in Fig. [6]). In Table [2]the rate of rock/ water r events are computed 
for the three different km 3 NT sites using several UHE neutrino fluxes as already 
considered in [HI [53] and described in [61]- [67] (see also Fig. [11]). For comparison, 
we also show in the last column the corresponding prediction for Earth-skimming v T 
at Auger-FD. As can be seen from Table [21 the relative enhancements due to matter 
effects on rock events can be as large as 30%, whereas the difference in the rates of water 
events for a fixed neutrino flux is mainly due to the different depth of the three sites. 

An interesting feature is the dependence of the event rate upon the shape of the 
NT detector for a fixed total volume of 1 km 3 , a property that might be relevant for 
the eventual design of the detector. Consider for example a km 3 NT placed at the 
NEMO site with the shape of a parallelepiped rather than a cube, where in particular 
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ANTARES 


NEMO 


NFRTDR 


Ang-pr-FD 141 1 


GZK-WB 


0.090/0.799 


0.107/0.929 


0.123/1.039 


0.074 


GZK-L 


0.099/1.076 


0.130/1.282 


0.157/1.465 


0.213 


GZK-H 


0.225/2.744 


0.313/3.280 


0.386/3.766 


0.560 


NH 


0.891/8.696 


1.102/10.19 


1.295/11.47 


1.245 


TD 


0.701/5.072 


0.817/5.799 


0.921/6.424 


0.548 



Table 2. Yearly rate of rock/water r events at the three km 3 NT sites for different UHE 
neutrino fluxes. GZK-H is for an initial proton flux oc 1/E, assuming that the EGRET 
flux is entirely due to it- photoproduction. GZK-L shows the neutrino flux when the 
associated photons contribute only up to 20% in the EGRET flux. GZK-WB stands 
for an initial proton flux oc 1/E 2 |61j-[65j. The other two neutrino fluxes correspond 
to more exotic UHECR models. NH represents the neutrino flux prediction in a model 
with new hadrons [66] , whereas TD is the neutrino flux for a topological defect model 
|67j . In the last column we report the corresponding prediction for Earth-skimming v T 
at Auger-FD. 

the E and W surfaces are enlarged by a factor 3 in the horizontal dimension, the N 
and S surfaces being reduced by the same factor, keeping the height of towers still of 
1 km. In this case the expected rate of rock events per year is enhanced by almost a 
factor 2, from 0.11 to 0.18 for the GZK-WB flux, while this enhancement could be even 
larger for neutrino fluxes with a larger high energy component. Moreover, the expected 
rate of water events increases as well by a factor of the order of 50%, from 0.93 up 
to 1.40 per year. Similar exercises can be also performed for the ANTARES and NESTOR 
sites. Of course, a further possibility which might favor UHE-r detection consists in 
increasing the effective volume of the detector keeping unchanged the 1 km height and 
the number of towers of photomultipliers but adopting a larger spacing. As an example, 
for a factor four larger volume with a doubling of the tower spacing both the rock and 
water r events would increase by almost a factor two, but obviously at the expense of 
the energy threshold and the quality of the event reconstruction for "low-energy" (TeV) 
neutrinos. For a detector aiming at the exploration of the range above the PeV, this is 
a less severe problem. 

The fact that the event rate depends upon the total surface of the detector is a 
peculiar feature of a NT, quite differently from what expected at the Auger observatory. 
Actually, in this case observed showers are generally initiated not very far from the 
detector compared to its dimensions so that the shape of the detector (i.e., the position 
on the border where the FDs are placed) is not as important as its "volume" (controlled 
by the area enclosed by the FDs). In fact, in order to produce a r emerging from the 
Earth with enough energy to generate a shower detectable by the Auger FDs, the energy 
of the neutrino should be larger than 1 EeV= 10 18 eV, taking into account the r energy 
loss in the rock. But the decay length of such a UHE r is /decay — 50 km x(£' r /EeV), 
to be compared with the dimensions of the Auger fiducial volume, ~ 50 x 60 x 10 km 3 . 
Conversely, a neutrino telescope can detect taus or muons which are produced very far 
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from the detector by a neutrino charged-current interaction, from distances comparable 
to the charged lepton range at that particular energy [31] . Indeed, the r range in water 
is of the order of several kilometers: from the value of f3 T = 0.71 x 10~ 6 cm 2 g" 1 we 
obtain an attenuation length l/(/3 T Q w ) — 15 km, while for muons (see next section) the 
range is approximately eight times smaller, of the order of 2 km. In other words, the 
effective volume of a NT of the kind discussed so far can be much larger than 1 km 3 , 
thus maximizing the detector area might greatly improve the detection rate. 

Of course, one should not forget that the design of a NT also depends strongly upon 
more detailed experimental considerations. Shapes which are not very compact or a 
detector with very sparse instrumentation have worse performances in the reconstruction 
of track properties as well as in signal-background separation, though this is mainly 
problematic at energies lower than 100 TeV, in the atmospheric neutrino energy range. 
In any case, our analysis suggests that the choice of the detector shape could be 
an important feature in orienting the target of a NT investigation towards either 
atmospheric or extra- atmospheric neutrino physics. In this respect, the possibility to 
modify this parameter quite easily for a NT water detector offers a great advantage with 
respect to an under-ice detector. 

A comment is in turn regarding the various approximations we used in the 
calculation. In particular we neglected tau regeneration effects. It is well known that 
these effects depend on the adopted incoming spectrum with a typical behavior in which 
the steeper the spectrum the less relevant the effect is [58| 134"] . For an E~ 2 spectrum 
the effect is almost negligible [58] . while for harder spectra like the ones we presently 
consider the effect can be of the order of 20% for taus coming from the nadir direction 
and of decreasing relevance for more horizontal events [34]. An estimate of the direction 
averaged effect gives then a correction less of 10%. 

A further approximation regards the stochastic nature of the tau interaction 
in matter that we approximated like a continuous energy loss process through the 
parametrization of Eq. ([51). At energies larger than 10 6 GeV, tau energy losses 
are affected by the large theoretical uncertainty on cross-section for photonuclear 
interaction, the leading mechanism at these energies (see [53] and [72], and references 
therein). Differently from muons, for taus the dominant source of uncertainty is not the 
stochastic vs. continuous nature of the energy loss, but the model-dependence of the 
photonuclear interaction itself; the stochastic nature of the losses is then sub-dominant 
with respect to the understanding of the process. A detailed discussion of the problem is 
given in Ref . [53] . The continuous approximation is then enough for the estimate of the 
mean rate values as given in the text, especially for a relative comparison of the sites. 
A more careful treatment would be needed if one is interested in a realistic estimate of 
the errors. 

A final comment deserves the dependence that matter effects could have not only 
on truly differences of amount of matter present in different directions, but also on 
the differences in the lepton interactions in water or rock due to the different chemical 
composition (A, Z), i.e. to differences between /3 r and (3 W . We study the issue performing 
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Figure 9. (Left) The effective apertures A T{ - r ' w) (E v ) for the NEMO site for a JA--nucleon 
cross section multiplied by a factor 0.5, 1 and 2 with respect to the standard result a. 
(Right) Ratios of the number of events in rock/water when the cross section is rescaled 
by the factor shown on the x axis, for several incoming UHE neutrino fluxes. 

detailed calculations of the lepton propagation as given in [73] . The calculations show 
that the (A,Z) dependence is of the order 10% for taus, almost constant at high energies 
(> 1 PeV), while of the order 20% for muons again almost constant at high energies. 
Given the model uncertainties in the tau losses and the level of approximation of 10% 
used throughout the paper we used the same value of (3 for rock and water losses for 
both muons and taus so that the differences seen at high energies in the apertures in 
Figj6]has to be all ascribed to a genuine matter effect. Note moreover that the matter 
effect is a feature increasing with energy, amounting for example in a factor even of four 
in difference between the W — E surfaces for NEMO at 10 11 GeV, while the percentage 
difference in (3 w /f3 r always remain of the order 10% through the whole energy range. 
The role of chemical composition then is at most subdominant. 

One of the main motivations for studying UHE neutrinos is that they provide 
a possibility to explore a range of energies for scattering processes which is still 
untested (maybe impossible to test) by particle accelerators. In this respect, measuring 
the neutrino-nucleon cross sections at high energies could have a large impact on 
constraining or discovering new physics beyond the standard model (see e.g. [301 EE])- 
While a measurement of the event energy spectrum cannot remove in general a 
degeneracy between the neutrino cross section and the incoming neutrino flux, a neutrino 
telescope could offer the interesting capability of disentangling these two factors because 
of the role of matter effects. Indeed, provided that enough statistics is collected and the 
detector has a good zenith angle resolution, the flux dependence can be subtracted off 
by measuring the ratio of the event rates coming from different directions [691 1701 171] . 
In the left panel of Fig. [9] we show how the NEMO effective apertures for r rock and 
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Figure 10. The effective apertures Aa (E v ) versus the muon neutrino energy for 
(left) rock events and (right) water events at the NEMO site. 



water events change if the neutrino-nucleon cross section is half or twice the standard 
model result, while in the right panel we display the ratio of water /rock r event rates 
for several adopted fluxes. We see that this ratio is quite sensitive to the value of 
the cross section. In particular, the number of rock events is essentially unaffected 
while the water event rate increases almost linearly with the cross section. Clearly, 
since the statistical error on the ratio would be dominated by the rare rock events, 
an experiment which aims at exploiting this effect should maximize the acceptance for 
almost horizontal events. We conclude by noticing that our results do not take into 
account any detailed experimental setup and up to this point the contribution is not 
yet considered. Nevertheless, since both the incoming neutrino flux and cross section 
on nucleons are expected to be flavor independent the possibility of determining both 
these quantities at a NT seems an interesting perspective. A more detailed analysis of 
this issue will be addressed elsewhere. 

4. The contribution: disentangling z/'s from r's 

In the previous Section we have discussed the rate of v T events implicitly assuming that 
a t lepton can be distinguished from a \x in a NT. However, given the experimental 
characteristics of the detector, this could be a difficult task and events should be also 
included in any realistic simulation. We address this issue in the present Section. 

Using the definition of the aperture in full analogy with Eq. ( 1121) and applying the 
same considerations of Section 2 to v^j /x, we have computed the \x apertures for water 
and rock events for the various surfaces of the NT, adopting the value (3^ = 0.58 x 10 -5 
cm 2 g _1 in the expression of the muon differential energy loss analogous to Eq. (jSJ) (as 
for the r, the term weighted by 7^ is negligible for the energy range of interest). The 
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Figure 11. The total NEMO sensitivities for (left) r and (right) [i events versus the 
neutrino energy, compared with the UHE neutrino fluxes considered in the paper. 

results are shown in Fig. [10J The main features as well as the role of matter effects are 
essentially unchanged for muons, the only difference coming from the muon contribution 
at lower energies because of the longer muon lifetime compared to that of taus. 

It is worth to discuss briefly the adopted parametrization for the muon energy losses. 
As discussed above, tau energy losses are affected by the large theoretical uncertainty on 
cross-section for photonuclear interaction. On the other hand, photonuclear interactions 
are less relevant for muon propagation, thus the theoretical uncertainty on the energy 
loss is correspondingly smaller. We then checked the validity of the approximation given 
in Section H] versus the detailed calculation given in [73] . We found that the accuracy 
is at the level of 15% over the whole energy range. The impact of this uncertainty 
on the expected event rate is as follows: a 15% increase of [3^ gives a few % decrease 
of the number of water events and a ~ 10 % decrease of the number of rock events. 
This uncertainty then does not affect the estimate of the number of events while a 
more careful treatment is required for a reliable forecast of the neutrino cross-section 
sensitivities at a Neutrino Telescope. 

In Fig. [TT]we summarize the r and \x results by showing the total sensitivity S ,M,T for 
NEMO defined as S^ T E U A^ T = 1 event year -1 (E u decade) -1 , with A^ T the total ji and 
r effective aperture, respectively; we also show the various neutrino fluxes considered 
through the paper. We see that in agreement with the results of Table El at least one 
event per year is expected even in the case of a GZK-WB flux, while larger rates are 
expected for higher fluxes (see also [71]). Notice that in the energy bin 10 8 — 10 10 GeV 
both fi and r contributions are comparable while muons are expected to dominate in 
the lower energy range, depending on the particular flux we consider. 

Concerning the possibility to distinguish between UHE taus and less energetic 
muons a comment is in turn. As we mentioned in the Introduction, the main difficulty 
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Figure 12. Total r, /j, and contained events apertures for the NEMO site. Note that 
the contained events are the same for all neutrino flavors. 



is that Cherenkov detectors like a NT do not measure the particle energy but rather 
the energy loss inside the detector volume and thus a high-energy r tracks can be 
misidentified with a muon track of lower energy. This is because the ratio of r to fi 
energy loss rate is given by Pr/Pn — 1/8. In principle, given that tau energy losses are 
dominated by photonuclear processes versus radiative interactions for muons, it would 
be possible to distinguish a muon-track from a tau-track from the different hadronic 
content along the track. However, Montecarlo simulations indicate that NT's are poorly 
sensitive to this signature [51] . 

As far as the contained events are concerned (i.e. where the charged lepton 
production happens inside the instrumented volume), the telescope has, instead, realistic 
chances of flavor-tagging. As long as the energy-loss or decay range of the particle to be 
detected is small compared with the detector size, the event rate depends basically on the 
fiducial volume. This is always the case for neutral current events, which if detectable 
produce a localized hadronic shower from the struck nucleon, and for charged current 
f e -events, since the electron rapidly loses its energy. On the contrary, for energies above 
the TeV scale, u^s charged current events produce muons which are detectable as tracks 
several km away from the production point. The case of v T charged current events is 
yet different: for energies < 10 7 GeV, the boosted decay range of the tau particle is 
negligible with respect to the detector size and depth, and the event rate is determined 
by the instrumented volume, like for z/ e -events. On the other hand, for the typical 
spacings between strings/photomultipliers considered by current designs, above the PeV 
scale the boosted r decay length is larger than the spacing and one starts resolving the 
tracks, while below this energy r produce showers which differ from z/ e -events only by 
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the hadronic content, which is however challenging to tag. 

We show in Fig. [12] the total apertures for r and fi events at the NEMO site together 
with the aperture for r contained events (the same would apply for \x contained events). 
Indeed, contained events represent a sub-leading although non-negligible fraction of the 
total number of events, always of the order of 10% for \x and even greater for low 
energy r ( < 10 7 GeV) due to the very short ( < 1 km) decay length at these energies. 
Moreover, the contained aperture depends only on the neutrino interaction probability 
so that it can be considered also a reliable estimate of the e-induced showering events, 
neglecting the effect of the Glashow v e e~ resonance at 6.3 PeV. It may also be possible 
to identify lollipop events in which a r with energy larger than PeV produces a long 
minimum-ionizing track that enters the detector and eventually ends in a huge burst 
as the r lepton decays into a final state with hadrons or an electron. In this case, 
the final burst would be a direct measurement of the r energy while the energy loss 
along the track would be smaller than for a muon of the same energy. Probably, the 
cleanest signature of a r particle would be the detections of a double-bang event [75] 
in which a v T interacts inside the detector and the produced r decays in shower again 
in the detector, but the probability of such an event is extremely small. In any case, 
all these possibilities suffer from lower statistics as they all require that the interactions 
(showering or production) occur inside the detector, with a reduction of the effective 
volume down to 1 km 3 compared with the several km 3 effective volume for /i and r 
events which go across the NT fiducial volume. 

In view of these considerations, we conclude that, at least for the bulk of the events, 
the most viable strategy is to combine both muon and tau contributions and construct 
spectra depending upon quantities which are directly observable. The simplest choice is 
to consider the energy loss rate inside the detector which amounts to measure the track 
length and the total deposited energy and the arrival direction. In Fig. [13] we show the 
contours of expected \x and r event rates in terms of the zenith angle of arrival directions 
and dE/dx ~ —/3Eg w . For relatively low energy losses (3Eq w < 10 5 GeV/km the whole 
dominant contribution comes from muons, which therefore can be easily disentangled, 
whereas in the high energy loss tails the event distributions are almost the same for 
both neutrino flavors and one is forced to use the total z/ M + v T events as the input of 
any analysis of the data. 

5. Conclusions 

Ultra-high energy neutrinos represent one of the main targets for several experiments 
which adopt a variety of detection techniques. Among these, the optical Cherenkov 
neutrino telescopes deployed under water or ice look for the tracks of charged leptons 
produced by the high energy neutrinos that reach the Earth. In this paper, we have 
presented a new study of the performance of a km 3 neutrino telescope to be located in 
any of the three sites in the Mediterranean sea proposed by the ANTARES, NEMO, and 
NESTOR collaborations. Our main goal is to compare the performances of the different 
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Figure 13. Contour plots of the number of rock events (red full lines) and water 
events (black dashed lines) at the NEMO site in the zenith angle-di?/ dx plane for r (left 
panel) and /i (right panel) assuming a GZK-WB neutrino flux. The contours enclose 
65, 95 and 99 % of the total number of events. 

sites, keeping apart detector-specific features, like partial detection efficiency or the 
architecture adopted for the towers of strings. We concentrated instead on the details of 
the under-water surface profile of each of the three sites, using the data from a Digital 
Elevation Map. By generating a realistic and statistically significant sample of v T jr 
and z/^/yU tracks crossing the fiducial volume of the km 3 neutrino telescope, we have 
calculated its effective aperture to UHE u T and neutrinos and the expected number 
of events for different UHE neutrino fluxes, for both cases where the neutrino/charged 
lepton track is crossing the rock (denoted as rock events) or the water only (denoted as 
water events). Our results can be summarized as follows: 

• The impact of the site geography (or matter effects) on observables such as the 
"rock fraction" of the total event rate or asymmetries in the event direction can be 
important, particularly at high energies. 

• Even for a fixed instrumented volume, these matter effects can be enhanced by a 
suitable choice of the geometry of the telescope, maximizing the lateral surface of 
the fiducial volume. 

• The continental crust provides an absolute orientation, hence the matter effects 
may provide a mean for calibrating the pointing capabilities of the detector, even 
when no point-source identification is possible, like for diffuse cosmic fluxes at UHE. 

• We analyzed briefly the dependence of the rock and water events from the neutrino 
fluxes and the neutrino-nucleon cross section. We found that the ratio of rock to 
water events may provide an additional way to disentangle the two unknowns, in 
addition e.g. to the well-known zenith angle dependence of the rates due to the 
screening effect of the average spherical Earth. Although a detailed analysis would 
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be needed, we stress that this may be important for constraining neutrino-nucleon 
cross section at UHE energies, otherwise unaccessible at the Lab. 

• While below the PeV scale the aperture for muon tracks is one order of magnitude 
larger than the aperture for contained events, above ~ 10 8 GeV the numbers of 
muon and tau track events are comparable. We have briefly addressed the problem 
of whether it is possible to distinguish /i's and r's in the detector. Apart for the 
sub-dominant fraction of contained events with specific signatures, we stressed that 
a realistic prescription at UHE is to sum the bulk of /i and r events, the natural 
variables for describing the events being the arrival direction and the energy loss 
rate in the fiducial volume. 

The main conclusion one can draw from our analysis is that the optimization for a 
telescope aiming at the E >PeV region is significantly different from one whose target 
is the E ~TeV range: in the first case, the search is basically background free and even 
a relatively poor angular and energy resolution may be acceptable. The crucial goal is 
to maximize the event rates, and the discrimination among models may be based on 
"counts" and a very rough directional and energy binning of the events. In this respect, 
one should maximize the instrumented volume — compatibly with the experimental 
requirements for a meaningful reconstruction of the event — and also carefully design 
the geometry to maximize the lateral surface, in order to exploit the matter effects 
provided by the underwater profile. On the other hand, at TeV energies angular and 
energy resolutions are crucial to improve the very signal to noise ratio, and may help 
identifying point-like sources. At the same time, in the TeV range the matter effects we 
have stressed on in this paper are less relevant, and should not influence significantly 
the choice of the site. 
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